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Energy transfer to acoustic phonons is the dominant low-temperature cooling channel of electrons 
in a crystal. For cold neutral graphene we find that the weak cooling power of its acoustical modes 
relative to the heat capacity of the system leads to a power law decay of the electronic temperature 
when far from equilibrium. For heavily doped graphene a high electronic temperature is shown 
to initially decrease linearly with time at a rate proportional to n 3//2 with n being the electronic 
density. We discuss the relative importance of optical and acoustic phonons to cooling. 
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Introduction — Energy exchange between the electrons in 
nanoscale electronic devices and their environment is a 
key issue in the design of electronic circuits and will play 
a role in any future graphene-based electronics [l|, 
The dominant electronic cooling mechanism in nearly any 
solid state environment is energy transfer to phonons. 
Energy relaxation in a graphene sheet is dominated by 
transfer to the acoustic and optical phonon modes of 
its two-dimensional honeycomb lattice and to the opti- 
cal phonon modes of its substrate [H, 0] • In this Letter 
we address electronic energy relaxation in graphene with 
a focus on the intrinsic cooling channel provided by the 
acoustical phonons. 

Optical measurements are a particularly useful probe 
of energy transfer between electrons and phonons and 
have been employed in the past in studies of electronic 
cooling in quantum wells Similar measurements were 
recently performed on epitaxial graphene samples @, 0|- 
In a typical measurement electrons are excited to high 
energies using an optical pulse. The relaxation process 
of the hot electrons is then monitored using differential 
transmission spectroscopy. Although transport measure- 
ments do in principle provide an alternate way of study- 
ing interactions between electrons and phonons the 
resistivity contribution from acoustic phonons in typi- 
cal graphene samples is much smaller than the elastic- 
scattering residual resistivity contribution Q. Transport 
is therefore relatively insensitive to the electron-phonon 
coupling strength. Even for suspended graphene samples 
in which transport is nearly ballistic quantum resistance 
dominates over the phonon induced resistance |9|. On the 
contrary, cooling of hot Dirac quasi-particles is entirely 
due to phonons. 

Guided by experiment 0, 0] we assume that e-e inter- 
actions thermalize the system throughout the relaxation 
process. Given the lattice temperature T L and the elec- 
tronic density n, the cooling process is then characterized 
by a single time dependent function, the electronic tem- 
perature T e (t). The time dependent chemical potential fi 
depends on T e and n and is determined by the conserva- 
tion of the number of particles. As the hot electrons equi- 
librate the electronic temperature decreases approaching 



its equilibrium value T L . 

A unique situation arises for graphene in the neu- 
tral regime, T e S> fi. In a typical semiconductor this 
non-degenerate regime is reached only at high temper- 
atures. In nearly neutral graphene, on the other hand, 
this regime is accessed at nearly all temperatures of inter- 
est. We find that when T L < T e < 180K the electronic 
temperature satisfies a power-law decay law, 
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Here To in the initial temperature of the electrons and 
D is the deformation potential measured in eV. Here- 
after we use H = 1 and measure all temperatures in meV. 
Transport measurements have been able to bound the 
value of D between lOeV and 50eV however, more pre- 
cise limits on this important parameter have remained 
elusive [9(. We propose that Eq.([T|) be used to obtain an 
accurate value of D from experiment. 

This slow low-T e cooling of neutral-regime Dirac quasi- 
particles differs markedly from the very fast cooling ob- 
served when T e js high. Recent experimental^, 0, E3] 



and theoretical 1 1] work has demonstrated that very hot 
electron plasmas (T e ~ several hundred meV) cool sig- 
nificantly after several pico-seconds. We find that low 
temperature cooling is slower by more than three orders 
of magnitude. 

Theory of Temperature Dynamics — The assumption of 
rapid thermalization implies that 



d t T e = Q/C 
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where C = dr e £ is the electronic heat capacity, Q = dt£ 
is the electronic cooling power and £ is the energy 
density of the system. Using the Boltzmann equation we 
find that 
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is the collision integral. Here a = v,c labels the valence 
and conduction bands whose energies are ±v g k with v g 
being graphene's band velocity. The occupation of each 
band is given by the time dependent Fermi distribution 
function /£ = /(e kQ , T e (t), p(t)) and 
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is the transition rate between state ka and state p/3. 
The energy exchanged with the phonon heat bath in the 
transition is e^p = — e P( g. In Eq.© N q = N(ui q ) is 
the Bose distribution function evaluated at the phonon 
energy u q . For acoustical phonons the transition ma- 
trix element is w q P = D 2 q 2 (l + s a /3 cos9)/Apui q . Here 
Sa/3 = 1 f° r intraband transitions and s a p = — 1 for the 
interband ones, 9 — 6^ — 6 p is the relative angle between 
the incoming and outgoing momenta and p is the mass 
density of graphene. For optical phonons « g 2 where 
g ss 2v g /a 2 y/2pAu!o with a = 1.42A and A being the area 
of the graphene sheet [l3j]. 

We first consider Q a the energy transfer to the acoustic 
phonon bath which is accurately described by the linear 
energy dispersion ui q = cq. It is instructive to separate 
the energy transfer into a loss due to spontaneous emis- 
sion, Q s P ont ^ and a gain due to induced transitions, Q md . 
Straightforward manipulations of (0J lead to 
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To evaluate Q a we exploit the large mismatch between 
v g and the sound velocity c and evaluate Q a to leading 
order in c/v g <C 1. In the limit c/v g — > the scattering is 
elastic, only intraband transitions are allowed, and there 
is no energy loss. For small c/v g interband scattering 
remains negligible and we can approximate |p| by |k| 
when performing the sum over q in Eq. ([?])■ In this way 
we find that to 0(c/v g ) 4 
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Similar steps yield Q s P ont and the total energy loss 

Q a = -^(T e -T L ) + (!-/«)]. (9) 
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As expected the net energy loss vanishes when the elec- 
tronic temperature reaches the lattice temperature. In- 
terestingly to leading order in c/v g the energy loss is in- 
dependent of the sound velocity. Due to the absence of 
inter-band transitions the cooling power from electrons in 
the conduction band and from holes in the valence band 
are simply additive. 

As is evident from Eqs.([8]) and ([9]) the energy gain 
due to the induced transitions is negligible compared to 
Qspont a £ m gh temperatures when T e ^> T L . However as 
the system cools the difference between Q md and | Q s P ont | 
decreases, vanishing in equilibrium. The vanishing of Q 
in equilibrium is assured by the detailed balance condi- 
tion reflected by the collision integral expression, Eq.(E]). 

The cooling power due to the intrinsic optical modes 
is easily estimated when the electron-phonon coupling is 
approximated by a constant g and phonon dispersion is 
neglected. Setting the phonon energy to luq it follows 
from Eqs.gHS]) that 
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Here N e and N L are the Bose distribution functions eval- 
uated at the temperatures T e and T L respectively, and the 
factor \x{x — 1)| originates from the electronic joint den- 
sity of states. For neutral graphene T w 1/6 when T e <C 
loq. The dominant optical phonon bath in suspended 
graphene is likely intrinsic. In unsuspended samples scat- 
tering by substrate phonons may introduce additional op- 
tical modes which are considerably less energetic |3[. We 
focus below on the energy loss due to acoustical phonons 
and comment on the role of the intrinsic optical phonon 
modes only at the end of this manuscript. 
Neutral limit — First we consider the neutral regime for 
which p(t) -C T e (t). If the equilibrium value of the 
chemical potential is significantly smaller than the lattice 
temperature the system remains in the neutral regime 
throughout the entire relaxation process, otherwise the 
system will eventually exit the neutral regime as it ap- 
proaches equilibrium. 

The energy exchanged in a typical transition is T e 
implying a momentum transfer of T e /v g . The typical 
phonon energy is then T e c/v g justifying the quasi-elastic 
approximation for the scattering by acoustical phonons 
for all values of the electronic temperature. This situ- 
ation is in marked contrast with the typical scenario in 
metals in which e-ph scattering becomes highly inelastic 
below the Bloch-Griineisen temperature. 

In the neutral limit Eq. ( [9]) for Q a can be further sim- 
plified by setting p, to zero in the integral to obtain a 
value proportional to T e 4 . Combining this result with 
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Eq.Q and noting that the energy density per spin and 
valley in neutral graphene is £ = 3C(3)Tg /2nv 2 where ( 
is the Riemann zeta function we find that 



d t T e = -7 T 2 (T e - T L ) 
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where 7 = 7vr 4 D 2 /540C(3)p^ = 1.18 • 10 3 L> 2 {meV 2 ■ 
sec) -1 with the deformation potential measured in eV's. 
Both the cooling power and the heat capacity decrease 
as T e approaches T L . The cooling rate slows because 
the former decrease is faster. Far from equilibrium 
when T e > T L Eq.(l2|) is solved by (J). The temper- 
ature decays as a power-law with a characteristic time of 
To = I/27TQ. The simple relationship between T e and 
the energy density £ combined with Eq. ( QJ yield 
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where £q is the initial energy density of the system. Like 
the electronic temperature, the energy density decays to 
equilibrium as a power-law when the system is far from 
equilibrium. 

Near equilibrium, when T e > T L , we linearize Eq. (ll2p 
with respect to T e — T L to find that the electronic temper- 
ature decays exponentially to its equilibrium value with a 
characteristic time given by t l = I/7T 2 = 848/Z? 2 Tj 2 /is. 

Our results for the neutral limit are valid to 0(p/T e ) 2 
since the particle hole symmetry of the system implies 
that both Q a and £ are even functions of pt/T e . 
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FIG. 1: (Color online) Equilibration of the electronic tem- 
perature. The evolution of the electronic temperature T e is 
plotted for electronic densities of (up to down) 0.02, 0.05, 0.1, 
0.2, 0.3 and 0.4 -10 12 [cm" 2 ]. The lattice temperature is lmeV 
and the deformation potential is assumed to be 20eV. The 
solid line corresponds to the equilibration of T c in a neutral 
system. 

Highly doped limit — We now turn to study the equili- 
bration of hot electrons in the doped regime for which 
p(t) ^ T e (t). In this regime we use the Sommerfeld ex- 
pansion to approximate Q a and C and obtain 
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where j d = 3D 2 el / 4n 2 pv* = 0.133L> 2 n 3/2 meV/nsec 
with n being measured in units of 10 12 cm~ 2 . Far from 
equilibrium the high electronic temperature initially de- 
creases linearly with time at an energy rate given by 7^. 
Near equilibrium T e approaches T L exponentially at a 
rate of jd/T L . 

The doped regime resembles the typical metallic case in 
that the quasi-elastic approximation breaks down at low 
temperatures when the electronic temperature is below 
the Blocb-Gritneisen temperature Tbg — 2ck w . There- 
fore our results for the doped regime describe the en- 
tire equilibration process for systems in which T L > T BG . 
However, for cold doped graphene our results are valid 
only when T e (t) > T BG . 

General solution — In the general case when p(t) ~ T e (t) 
the time evolution of the electronic temperature and 
chemical potential follow from two coupled differential 
equations. The first equation 
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is obtained from the number equation dtn = and 
expresses the conservation of the number of particles 
throughout the relaxation process. Here 
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The second differential equation 
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FIG. 2: (Color online) Equilibration of the chemical potential. 
The evolution of the chemical potential p is scaled with the 
Fermi energy and plotted for electronic densities of (from bot- 
tom to top) 0.02, 0.05, 0.1, 0.2, 0.3 and 0.4 -10 12 [cm -2 ]. The 
lattice temperature is lmeV and the deformation potential is 
assumed to be 20eV. 

We have solved the coupled differential equations 
(|15ll7p numerically for various values of densities. The 
electronic temperature as a function of time is plotted 
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in FigED for T L = lmeV, T = 50meV and D = 20eV. 
Clearly the equilibration process is faster for doped sys- 
tems. 

The corresponding evolution of the chemical potential 
for the different electronic densities is plotted in Figf2] 
The conservation of particle number throughout the equi- 
libration process forces \fi\ to increase with time. 




FIG. 3: The energy loss due to optical phonons divided by 
the energy loss do to acoustical phonons for neutral suspended 
graphene is plotted vs. the electronic temperature for lattice 
temperatures of (bottom to top) 1, 10, 20, and 30 meV. 

Discussion — The small low-temperature cooling power 
of the neutral graphene electronic system is due to both 
the small joint density-of-states for electronic transitions 
and the small energy of acoustic phonons at typical tran- 
sition momenta. From dimension analysis of Q and £ and 
Eq.© we find that the instantaneous energy decay rate 
a = \d t ln(£)| = Q/\£\ of a d-dimensional gapless neutral 
system at zero lattice temperature with conduction and 
valence band dispersions oc k s satisfies 



work on the cooling of bilayer graphene. The evolution of 
the electronic temperature of a gapless 2D system with a 
parabolic dispersion is given by d t T e — — jy/T^ (T e — T L ). 
However we do not expect this relationship to apply pre- 
cisely in bilayer graphene since the momentum depen- 
dence of its energy spectrum can not be described by a 
single power 
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Although the main focus on this work has been energy 
loss due to acoustic phonons we emphasize that the en- 
ergetic optical phonons will play a dominant role in the 
cooling of graphene at sufficiently high temperatures. To 
estimate the regime for which acoustic phonons domi- 
nate cooling in graphene we plot Q a / Q a as a function 
of T e . We use a simple model in which ujq — 196meV 
for both longitidinal and transverse optical branches. In 
Figj3] we plot Q / Q a for a neutral graphene sheet for 
various values of T L . Surprisingly the different functional 
dependencies of Q and Q a on the lattice and the elec- 
tronic temperatures lead to a non monotonic dependence 
of Qo/ Q a on T e and to the dominance of Q near equilib- 
rium at moderate lattice temperatures. As evident from 
Figf4] the relative cooling power of the acoustic phonons 
increases with doping, increasing the maximum temper- 
ature at which they are dominant. 

The inefficient cooling of graphene by acoustic 
phonons has immediate consequences for its non-linear 
transport [H, . 
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Exponential decay of T e and £ occurs when a is constant, 
i.e. when s = d+1. In systems like single layer or bilayer 
graphene with s < d+1, the cooling rate will decrease as 
a power law of the instantaneous temperature when the 
system is far from equilibrium. 
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FIG. 4: The analogue of figure [3] for a doped system with 
n = 10 13 cm" 2 . 



It is interesting to consider the consequences of this 
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